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Chapter 1 



Macroscopic phase coherence and 
Josephson effect 

In this Chapter we consider systems, that at first glance may seem very different: (i) 
superconductors, which at temperatures larger than the critical temperature T c (T c in 
conventional superconductors is of order of a few K) turn metallic and are described 
by the Fermi Liquid theory, (ii) superfluid liquid He: bosonic 4 He, T c = 2.17 K and 
fermionic 3 He, T c = 0.0025-fT; and (iii) Bose-Einstein condensates of cold alkali atoms 
(with T c of the order of nano Kelvin). All these systems have one important property 
in common - at low temperatures they possess macroscopic coherence. 

In previous Chapters it was discussed that a superconducting state is a state with 
a broken U(l) -symmetry and is thus characterized by a complex order parameter 

^ = tf( r ) = |^( r )|e^ (r) (1.1) 

which represents a macroscopic wave function of a superconductor. The presence of the 
order parameter means that at any given moment the phase difference of ^/-functions 
between any two macroscopically separated points in the superconductor is fixed, so 
that the whole sample acquires macroscopic phase coherence, or in other words, long- 
range order develops. One can describe a condensate of cold bosons or a superfluid He 
system in a similar way, since these systems are also phase coherent. 

Phase coherence leads to a number of specific quantum effects. For example, in 
superconductors it causes the quantization of magnetic flux first considered by London. 
One of the most celebrated manifestations of the phase coherence property is however 
the Josephson effect which is the subject of the current Chapter. 
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Chapter 2 

Josephson Effect in superconductors 



In 1962 Brian Josephson predicted a curious effect occurring in a system of two weakly- 
linked superconductors [T]. He demonstrated that a direct current can flow between 
two superconductors coupled via an insulating thin layer although no external voltage 
is applied. Furthermore, he showed that an external voltage would give rise to a rapidly 
oscillating current. Josephson's theory of a current induced between two superconduc- 
tors was rather fast confirmed experimentally [2]. 

The Josephson effect can be understood by the following simple considerations. 
First of all it is important to observe that a gradient of the phase gives rise to a current 

, = -^-V* (2.1) 

where N s is the number of superconducting electrons. Imagine now that instead of a 
uniform superconductor we deal with a superconductor with impurities: point-like, ran- 
domly distributed non-magnetic impurities. As we know, for an s-wave superconductor 
impurities do not affect the critical temperature, they do not destroy phase coherence 
in a superconductor (in accordance with Anderson's theorem [3]). As a consequence 
supercurrents (12.11) can flow through the system without a problem. This situation does 
not change if instead of impurities distributed in the bulk we have impurities distributed 
solely in a plane inside a superconductor (Fig J2.1( a)). 

Supercurrents proportional to the the fixed phase gradient can flow through the 
plane, even though the mean spacing between the impurities is smaller than the coher- 
ence length of a superconductor. This intuitively clear picture can be generalized to 
the case of a somewhat more complicated system: two superconductors separated by a 
thin insulating layer. This nontrivial generalization was realized by Josephson and led 
him to the prediction of two effects named after him: a.c. and d.c. Josephson effects. 

These Josephson effects are pure quantum phenomena, because electrons travel from 
one superconductor to the other by means of quantum mechanical tunneling through 
the barrier separating the two superconducting systems. The presence of a barrier, or 
inhomogeneity leads to the fact that the phase has a jump at the barrier (Fig. 12.1( c)). 
The supercurrent through the barrier is then driven by the phase difference </> = </> 2 — </>i . 
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Figure 2.1: Superconductor with an "impurity plane" inside (a). Supercurrent will per- 
sist even if the plane is replaced by a thin insulating layer (b). The phase jump at the 
barrier separating two superconductors (c). 



Thus, a supercurrent can flow between two superconductors provided they are separated 
by a sufficiently thin insulating layer. This effect is referred to as first or stationary, 
or d.c. Josephson effect. In this case the potential difference through the barrier is 
equal to zero. Note, that while the superconducting coherence length is of the order of 
£ ~ 10 4 A, the thickness of the insulting layer should be of the order of 10 — 20 A, i.e. 
thousands of times smaller than £. 

When a finite voltage bias V is applied to an SIS junction, the second, non-stationary 
or a.c. Josephson effect can be observed. In this case the current will be oscillating 
between the two superconductors with a frequency u proportional to the applied bias 

hw = 2eV. (2.2) 

Since in a superconductor electrons are bound into Cooper pairs, these pairs par- 
ticipate in the tunneling across a barrier between two superconductors. The energy 
2eV is then just the difference in the energy of a Cooper pair in passing from one 
superconductor to the other. 

One can estimate a Josephson current from simple electrodynamic considerations 
[4]. First of all we observe, that the appearance of a current in the system is related 
to the excess energy AE = E(<j>) — E(<fi = 0) associated with the junction between the 
two superconductors. It is clear that AE should be proportional to the product of two 
superconducting gaps AiA 2 , because if one of superconductors is absent, the excess 
energy vanishes. Apart from that the excess energy should be real, one can therefore 
conjecture its simplest possible form as follows 

AE = aj dxdy (j A X A 2 | - ~(AiA£ + c.c.) J = 2a J dxdy\ A 1 A 2 |(1 - cos </>), (2.3) 

where a is a phenomenological constant describing the coupling between the two su- 
perconductors in the x — y junction plane. 
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One can derive now from the gauge invariance principle that the supercurrent is 
proportional to (sin0). The gauge invariance requires the replacement 



he 



(2.4) 



where A is the vector potential. For simplicity we can chose the vector potential 
perpendicular to the x — y plane, so that integration along that axes of the r.h.s. of 
M gives 



7- ' 

he 



(2.5) 



where 4>i is the phase of the "left" superconductor, and 2 is the phase of the "right" 
superconductor. The excess energy becomes 



AE = 2a J dxdy\A x A^ 



cos 



2e f R , 
— I dzA, 

he 



(2.6) 



Variation of this energy with respect to the potential gives 



S(AE) = ~j~ a 



- |^ I 

he 



SA Z . 



(2.7) 



From electrodynamics we know that 



5E 



j 5 A dV, 



(2.8) 



and we get for the supercurrent 

2e 



j 



h 



a|A 1 A 2 | sin 0i - 2 - 



2e 
He 



R 



dzA, 



(2.9) 



For zero vector potential we obtain the famous expression for the Josephson current 

j = Jc sm{<f) 1 -<f) 2 ). (2.10) 

The current vanishes for 4>i ~ 4>2 = 0. The so-called critical current j c should be 
calculated microscopically (see Section EJ). 

The basic result (12.101) can be obtained in a different way. The following derivation 
in terms of a two-level system is due to Feynman [5]. As we discussed already, each 
superconductor can be considered as a macroscopic quantum state described by a wave- 
function \lZi 2 (11.11) . Since the coupling between the superconductors is very weak, the 
state vector describing the coupled system can be written in a simple form 



|$) = * 1 |1) + * 2 |2). 



(2.11) 
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The density of superconducting electrons in the left (right) superconductor, described 
by the state |1) (|2)) is defined as 

n t = \%\ 2 = (2.12) 

i = 1,2. The Schrodinger equation of motion for the state vector (12.111) reads 

ihd t \^) = H\V), (2.13) 

with the Hamiltonian 

H = H 1 + H 2 + H mt . (2.14) 

Here Hi = Ei\i)(i\, (i = 1,2) and the coupling between the superconductors can be 
written in analogy with ( 12.31) 

^ = -f(|l)(2| + |2)(l|). (2.15) 

Projection of the Eq. (12.131) on the two states gives the equations of motion for two 
weakly coupled superconductors 





= E^ x - 


a 

2* 2 < 




= E 2 m 2 - 


2 



(2.16) 

Remembering that can be expressed in terms of superconducting densities (see 

% = vW*, (2.17) 
we can derive the final equations in terms of the densities and phases 

a 

d t ni = -yjnin 2 sin(0i - <p 2 ) = -d t n 2 , (2.18) 

dt{<h-<f>l) = \{E, - E 2 ) + €03(0! - <P 2 ). (2.19) 

n In y/nin 2 

For equal densities rii = n 2 = n we get 

d t ((p2-(pi) = ^(E 1 -E 2 ). (2.20) 

The pair current density is given by 

j = dux = -dn 2 = j c sin(0i - <p 2 ). (2.21) 

with j c = an/h for equal densities. One should note, that the densities ri\ and n 2 are 
considered to be constant (we will see that this is not the case in a Bose Josephson junc- 
tion), their time derivative is however not constant due to the presence of the external 
current source which continuously replaces the pairs tunneling across the barrier. 
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The presence of the potential difference V is easily taken into account in our equa- 
tions. In two isolated superconductors the energy terms are given by the chemical 
potentials Ei = 2/ij {i = 1,2). A d.c. potential difference will shift these chemical 
potentials by eV, so that E\ — E 2 = 2eV, and Eq. (12.191) becomes 

2eV 

dt{<h - 0i) = (2-22) 

The two equations (12.2ip and (12.221) constitute two main relations of the Josephson 
effect, which we discussed at the beginning of this Section. For V = the phase 
difference is constant, so that a finite current density with a maximum value j c can flow 
through the barrier with zero voltage drop across the junction. This is the essence of 
the d.c. Josephson effect. With a finite potential difference V applied to the junction 
there appears an alternating current 

j = Jc sm(4> + jVt) (2.23) 

with a frequency (12.21) . This corresponds to the a.c. Josephson effect. 

In the following we derive the microscopic expression for the critical Josephson 
current. 



Chapter 3 



Microscopic derivation of a critical 
superconducting current 

The microscopic approach for the calculation of critical current was suggested by An- 
derson [Sj, and Ambegaokar and Baratoff [7]. Their method is based on the so-called 
tunneling Hamiltonian. In this approach the details of the interface are not taken into 
account and instead two weakly coupled superconductors described by Hamiltonians 
Hi and H 2 in the absence of tunneling are considered, whose coupling in the first order 
perturbation theory is described by the tunneling term in the Hamiltonian 

H = H 1 + H 2 + H T . (3.1) 

Ht has a simple form 

H t = Y1 T pq a pA- + T p A«p-> (3-2) 

pqcr 

where a are the fermionic operators of the "left" superconductor, and b are the fermionic 
operators of the "right" superconductor, a is a spin index, p and q are the momenta 
of electrons. Due to the time-reversal invariance of the Schrodinger equation (t — > — t, 
\jz _> the matrix elements T pq have the property T pq = T^ p The Hamiltonian 
(13.21) conserves the number of particles in the system N\ + N 2 , where 

Ni = Yl a P- a P- N * = Y1 6 pA- (3-3) 

per per 

The current is related to the change of the number of particles with time and is 
therefore by definition 

I = e(N 1 ) = -e(N 2 ). (3.4) 
The equation of motion for the operator JVj reads 

ihN i = [N h H T ), (3.5) 
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where we took into account that the operator AT* commutes with Hi and H 2 (i — 1,2). 
With (|3.5p the expression for the current (13.41) becomes 



1 = ~J [ ^2 T Pi( a U b w) ~ ^AaOpo) J ■ ( 3 - 6 ) 

One can proceed with the derivation of the Josephson current in several ways, see for 
instance pjj E]. Here we suggest a rather straight-forward derivation based on the 
nonequilibrium Keldysh technique [TOl [IT]. We apply then our general results to a 
problem of a stationary Josephson current between two superconductors. 

ilra t 




Figure 3.1: The Keldysh path on the time-plane along which the nonequilibrium Green's 
function is defined. 

One can introduce the so-called Keldysh Green's function [ID] which is defined in 
the following way 

G* CT (MO = 2i(at CT (OMO>- (3-7) 
The Keldysh Green's function ( 13.71 ) is a part of a general path-ordered Green's function 

CW(M') = -i(T c M*X(0>, (3-8) 

which is defined on a so-called Keldysh contour C = G\ + C% showed in Fig. 13.11 Tc is 
a time-ordering operator along this contour. It is convenient to separate the Keldysh 
contour on the upper (C\) and lower (C2) contours and to present the function ( 13.80 in 
a matrix form 

G=(r U rO' ( 3 - 9 ) 
\ ^21 ^22 / 

where indices 1 and 2 refer to the upper or lower Keldysh contour respectively. By 
applying the usual rotation operator 

h =^(-i0 (3 - io) 
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to the matrix (13.91 ): G — > R 1 GR we get the nonequilibrium Green's function in the 
Larkin-Ovchinnikov representation |12| 

G=( „ G ,j. (3.11) 

The relation between Gn, G12, G21, G22 and retarded G R , advance G A and Keldysh 
G K Green's functions can be trivially derived from (13.91) and (13.101) . For details see 

The Keldysh Green's function (13.7ft is useful, because we can immediately express 
the current (13.61) in terms of such functions 

7 = E (^SU*) - ^pvw) • (3-i2) 

pqcr 

One can readily see that [G pcia (t)}* = — G^ pa (t), so that the expression for the current 
becomes even simpler 

(3.13) 



h 



y^^pq^qpa^) 

We need thus to calculate the Keldysh Green's function G^ (t). For simplicity we 
proceed with our calculations in the first order of perturbation theory. We also consider 
only G pq | function, because G pq | can be derived analogously. 

The Green's function of the system described by the Hamiltonian (13.11 ) in the first 
order of perturbation theory reads 

G&M = -*( T c (l -ij drH T {r)\ b^{t)a^(t')). (3.14) 

According to the BCS theory only electrons with opposite momenta and spins are 
allowed to pair, hence we get the expression 

G« T (M') =T_ p ^(T c [ dr6 qT (06- qi (T)al pi (r)a pT (t')>. (3.15) 

Jc 

We assume, that in the simplest approximation the critical supercurrent is carried by 
Cooper pairs, and by applying Wick's theorem [13] to ( 13.151) we get 

Gpq T (*>*') = -r- P ,- q / dT^ u (q,t-r)T^(p,T-t'). (3.16) 

Jc 



Here we introduced anomalous Gor'kov functions for a superconductor [1 3j 

£=u(q,f-0 = -*<7bM^-q#)> 

T n (p,t-t') = -i{T c al pl (t)a^(t')). (3.17) 
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These Green's function do not depend on the sign of momentum p, but the order of 
spin indices does matter (replacing |j with JT will give an extra minus sign). 

In a lengthy but straightforward calculation [TT] one can extract the Keldysh part 
of the matrix Green's function (I3.16P 

/oo 
dr[^(q, t - r)Ff T (p, r-t') + ** (q> t - r)T*(p, r - t% 
-oo 

(3.18) 

The Fourier transformation of this expression gives 

[G { Z(")] K = -^[^(q,w)^f T (p,«) + ^(q,w)^(p,«)]. 
The current then becomes 



(3.19) 



2e 



Re 



pq 



(3.20) 



the factor of "2" appears because the contribution from G^^ is equivalent to the con- 
tribution from Gp ^. 

This expression simplifies greatly in equilibrium, in which case the Keldysh Green's 
function can be expressed as 



G*( W ) = (G»-G»)(l-2/( w )) ) 



(3.21) 



where f(u) = [e^^ + 1) 1 is the Fermi distribution function. Hence we get 



2e 

/ = —Re 
n 



pq 



T^tanh [^j (^S(q,w)^ft(p,w) - ^fj(q, ^^(p, u] 



(3.22) 

Substituting the explicit expressions for the retarded and the advanced Gor'kov 
functions [T3l we obtain 



2e 



pq 



du ( u \ 

tan h I I 



2tt 



-A? 



V277 V(w + i5) 2 -e\- |A q | 2 (a, + i5) 2 - e 2 p 



Ar 



A r 



{u - id) 2 -e\- |A q | 2 (u - id) 2 -el - |A, 



(3.23) 



We assume that superconducting gaps in the left and right leads are momentum- 
independent A q = Ai and A p = A 2 , and their product gives A1A2 = | Ai| | A 2 |e 8,i , 
where is the phase difference between two superconductors. We also take into ac- 
count that the integrand is purely imaginary, so that the current becomes 



2e 



lAxHAaKisin^r^lTpql 2 / 

L P q J ° 



du 



1 



1 



2vr \(u + i6)' 



IZS.il 2 (u + iS) 2 -el 



(lu - i5) 2 - el- \A!\ 2 (cu-iSy 



p 

12 



A, 



tanh 



(-) 

\2TJ 



■|A 2 | 2 
(3.24) 



We thus derived microscopically that the current between two superconductors is pro- 
portional to sin (p 

/ = / c sin0, (3.25) 

where I c is the critical supercurrent 

2«e, A ,, A , „ f 00 duo ( 1 1 



pq 



27r\(uj + i5) 2 - el - |Ai| 2 (cj + 25) 2 -e p - |A 2 



1 , / U) 



tanh — . (3.26) 



(u - id) 2 -e\- |Ai| 2 {uj - id) 2 - e p - \A 2 \ 2 J \2T J ' 

Using standard rules of contour integration we can replace the integral over uj by a sum 
over discrete Matsubara frequencies uj n = (2n + 1)ttT [9] 

4 °° 1 1 

I c = -^|Ai||A 2 |T V|T pq | 2 V — - 5 r-r-^. (3.27) 

h ii n f^^ + e q +|Ai| 2 uj 2 n + e 2 p +\A 2 \ 2 

The summation over p and q can be replaced by an integral, which can be taken 

a 2 °° 1 1 

h= ^|A 1 ||A 2 |A^ 1 (0)^ 2 (0) \T \ 2 T V , (3.28) 



where iVi(O) and iV 2 (0) are the densities of states at the Fermi energy in the normal 
state of left and right lead correspondingly. When the two gaps are equal to each other 
Ai = A 2 = A, this expression takes a simple form 

I c = ^A iV\(0)iV 2 (0) |T | 2 tanh (J^j . ( 3 .29) 

One usually introduces the so-called resistance of the tunneling junction in the normal 

state Rn 

i- = 47^(0)^(0) |T | 2 , (3.30) 

tin 

so that 

ir A ( A \ 
7 cJ R n = — tanh(^— J. (3.31) 

In Fig. 13.21 we depict the temperature behavior of the Josephson critical current I c , 
normalized by 7iT c /2eR n . In order to obtain this dependence we had to solve the 
standard BCS gap equation, and we present the temperature dependence of the gap in 
the same Fig. One can see, that the critical current is monotonously decreasing with 
temperature rather similar to the gap behavior. Near T c the current is proportional to 
A 2 and is therefore linear in (T c — T) . 

At zero temperature the critical supercurrent for superconductors with different 
gaps is 

2A 1 (0)A 2 (0) /|A!(0)-A 2 



4 ei?[Ai(0) + A 2 (0)]^ \ Ai(0) + A 2 (0) ' ' (3 ' 32) 

13 



where K(x) is a complete elliptic integral of the first kind. In case of Ai(0) 
A we get a simple expression 



= A 2 (0) = 
(3.33) 



I 




Figure 3.2: Josephson critical current defined by Eq. (13.311) in the units of irT c /2eR r 
(solid line). The dashed line shows the temperature dependence of the BCS gap A/T c . 
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Chapter 4 



Josephson effect in Bose-Einstein 
condensates 




Figure 4.1: Bose Josephson junction: two Bose-Einstein condensates are trapped in an 
external double-well potential V ext . 

Condensates of cold alkali atoms [HJ [15] provide a unique opportunity to realize 
and to control Josephson effect in a weakly interacting bosonic system. It was first 
predicted in [H]. A weak link between two condensates can be realized in a double- well 
external potential V ext (Fig. 14. ID . An interacting system of bosons confined in such a 
potential is described by a general Hamiltonian 



H 



dr \j/ f (r,t) 



2m 



V 2 + V ext (r 



*(r,i) 



dr / dr'&(r,t)&(r',t)V(r-r')ty(r',t)ty(r,t). 



(4.1) 



Here ^(r, t) is the bosonic field operator and V{r — r') is a two-particle interaction. 

At low temperatures an experimentally realized gas of bosons is very dilute and par- 
ticles are weakly interacting, one can therefore introduce a contact interaction between 
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the particles 

V(r-r') = g6(r-T r ), (4.2) 

where g = Airh 2 a s /m with a s being the s— wave scattering length. According to the 
Bogoliubov approximation one can consider the bosonic field operator as a sum of a 
classical field (condensate wave function, representing the condensate order parameter) 
and excitations 

$(r,t) = *(r,t) + <J*(r,t). (4.3) 

In the mean-field description we can neglect the excitations due the smallness of the 
interaction term, so that our Hamiltonian becomes essentially classical 



H = dr &(r,t) 



h 2 

- — V 2 + V ext (r) 
2m 



tf(r,t) + | J dr^(r,t)^(r,t)^(r,t)^(r,t). 



(4.4) 

For the Josephson effect to occur only a small overlap of the condensate wave-functions 
is sufficient, and we can assume that the condensate wave-function \I/(r, t) is given by 
the sum of the order parameters for each well |17| 



(r, t) = ^x(r, t) + tf 2 (r, t) = <^i(r)^i(t) + ^ 2 {v)^ 2 {t). (4.5) 
Here y?i(r) and tp 2 (r) are the ground state solutions for isolated traps [TZI [THJ, and 

Vi(t) = Ni(t)e^ (4.6) 

is the complex condensate order parameter with Ni being the number of particles in the 
i-th well, and fa is the phase of the condensate in the same well. With these notations 
the Hamilton function (|4.4I) takes the form 

H = EkN! + E 2 N 2 + ^jN 2 + yiV 2 2 + 2J y /N 1 N 2 cos((l>), (4.7) 

where <f> is the phase difference between the wells, 

Ui = g [ rfr|^| 4 , (4.8) 



Ei = j dx(^\V^v)\ 2 + ^{v)VU^)\ (4-9) 



and J is the Josephson coupling 



J = - J dr ^A_vy>i(r)Vp 2 (r) + ^i(r)^ 2 (r)K*t(r)^ . (4.10) 

One can reexpress the Hamiltonian (14.71) in terms of a particle imbalance n 

N 1 -N 2 , 

71 = -r± rf 4.11 

^+N 2 ' 
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so that the effective (dimensionless) Josephson Hamiltonian reads 



A 



H= -n 2 - y/i - n 2 cos + AE n. (4.12) 

We see that in this case only two parameters determine the behavior of the system: the 
effective interaction 

A = (E/i + ?7 2 )(iVi + N 2 )/{AJ) (4.13) 
and the effective chemical potential difference 

AE= E.-E 2 {Ul -U 2 )(N 1 + N 2 ) 

2J 4J y ' ' 

As the particle number operator n and phase operator are canonically conjugated 
variables, in the classical case one can identify a corresponding Poisson bracket with 
their commutator. We can then derive the corresponding equations of motion for the 
particle imbalance and phase difference 

dH ■ dH .„ . 

n = ~W (415) 

As a result we get 

n = —\J 1 — n 2 sin 0, (4.16) 

n 

= An + AE + cos0. (4.17) 

V 1 — n 2 

These equations can be solved exactly in terms of the elliptic functions [20]. In Fig. 
14.21 we show the numerical solutions of the equations ( 14.171) for AE = 0. One can 
observe a qualitative change in the oscillations after A exceeds a certain "crossover" 
value A c . For initial conditions as in Fig. 14.21 A,. ~ 10. For A < A c the time-average of 
the particle imbalance is zero: (n) t = 0. For larger values of A the particle imbalance 
oscillates around a finite value (in case of Fig. 14.21 (n) t = 0.4). It means that on 
average the number of particles in one well is larger than the number of particles in 
the other well. This curious quantum phenomenon was termed macroscopic quantum 
self-trapping (MST) [17]. This behavior can be also seen in a phase portrait in Fig J4.3l 
showing the constant energy trajectories for different A. The "running" trajectories (c) 
and (d) correspond to MST. 

The occurrence of the MST phenomenon is readily understood if one remembers 
that the canonical Josephson Hamiltonian (14.120 in the small n limit 

H = ^n 2 - cos + AE n (4.18) 

can be mapped onto a pendulum Hamiltonian with tilt angle 0, dimensionless angular 
momentum p<f> = n, inverse mass A and applied torque AE [15]. For small n the bosonic 
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10 15 

2Jt 



Figure 4.2: Temporal oscillations of the particle imbalance n and phase difference 
A = 5 (a), A = 9 (b), and A = 11 (c). Initially n(0) = 0.6 and 0(0) = 0. 



or 
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-0.2 








-0.4 






-0.6 







Figure 4.3: Phase space plot: (a), (b) and (c) for same parameters as in Fig \4-2\ (d) is 
for A = 11, n(0) = -0.6, 0(0) = 2tt. 



junction supports small-amplitude Josephson "plasma" oscillations with the frequency 
uj = y/X/(2J) Fig. 14.4( a). A rotation of the pendulum in Fig. 14.4( b) corresponds to 
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Figure 4.4: Analogy with pendulum oscillations: MS T state corresponds to the case (b), 
while standard Josephson oscillations of a bosonic junction to the case (a). 

the MST-state with a running phase. It is also easy to understand the "wiggles" in the 
phase dynamics in Fig. 14.21 (c). as the pendulum always slows down at its highest point. 

The physical behavior of the Bose Josephson junction described by the Hamiltonian 
(14.121) is more complicated due to the additional factor yl — n 2 . In the language of 
pendulum analogy it means that the length of the pendulum is not rigid anymore, but 
varies with time. This leads to additional fixed points in comparison with the previous, 
simple example (for a detailed description see [20l fT9] ). 

The value of A c which determines the crossover to the MST state is defined by the 
condition 

Ho = H(n(0), 0(0)) = ^n(O) 2 - ^1 - n(0) 2 cos(0(O)) > 1, (4.19) 

so that 



l + yl-n(oy»oosMO)) 

It means that A c can be relatively easy controlled in an experiment by varying the 
initial conditions for two condensates. This property was used by experimentalists and 
the predicted in [T7] behavior of the Bose Josephson junction was successfully observed 
experimentally [21] . Josephson oscillations of the particle imbalance with (n) t = 
were observed for n(0) ~ 0.28(6) and 0(0) ~ 0. The MST regime was achieved with 
n(0) ~ 0.62(6) and 0(0) ~ 0, the parameter A in both cases is estimated to have a fixed 
value of 15(3). 

The analogue of a.c. and d.c. Josephson effect discussed for superconductors in 
Section [2] has been recently observed in a Bose Josephson junction |22j. 
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Chapter 5 

Josephson effect in superfluid He 



Since superfluid Helium possesses phase rigidity, one would expect the Josephson effect 
to occur between two weakly coupled Helium systems. Although liquid Helium was 
discovered more than seventy years ago (in 1937 by P. Kapitsa, J. F. Allen and D. 
Misener), it took a long time before the Josephson effect has been finally observed: 
in 1997 in superfluid 3 He [23] and in 2001 in 4 He |24| . The main obstacle for the 
observation of the effect is a very small coherence (healing) length of Helium: 50 nm 
for 3 He and even smaller, of the order of 0.1 nm for 4 He. It took thus almost 60 years 
to overcome two main technical difficulties: (i) the creation of the weak link itself - a 
structure with small apertures with dimensions of the scale of the coherence length, (ii) 
measurement of tiny mass currents which would flow through such a structure due to 
the Josephson effect. 

Note, that in case of liquid Helium, one can not apply an external electromagnetic 
potential difference to the system, as in the case of superconductors, neither can one 
modify the trapping potential to simulate this effect as in Bose-Einstein condensates 
of cold atoms. For Helium the role of external potential V is played by pressure P, so 
that both Josephson relations (12.211) and (12.221) remain the same with chemical potential 
difference proportional to pressure: 

APm , . 

Ap = ——. (5.1) 

Here m is the mass of either the 4 He atom, or twice the 3 He atomic mass ( 3 He is a 
fermionic system and its superfluidity is induced by coupled fermions), p is the liquid 
density. Applied pressure difference will induce therefore an oscillating mass supercur- 
rent with the frequency uoj = APm/ ph. 

In the experiment [23] two 3 He systems are separated by a membrane with numer- 
ous apertures only 100 nm in diameter. The healing length at the given experimental 
temperature was slightly exceeding the aperture diameter. The great number of aper- 
tures (more than 4000) served to coherently increase the monitored supercurrent, which 
was otherwise too tiny to be resolve in the measurement. Another soft membrane con- 
trolled by an applied bias was used in order to create an external pressure difference. 
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Any displacements of the membrane due to the supercurrent were monitored. Finally, 
the signals of the supercurrent frequency obtained in this way were amplified and con- 
nected to audio head-phones, and the listener could literally hear the effect of coherent 
quantum oscillations between weakly coupled superfluids. It sounded like a whistle 
smoothly drifting from high to low frequency while the pressure relaxed to its zero 
value. The dependence of the supercurrent frequency on AP has been found to be 
perfectly linear |23j. 

In 4 He the regime of ordinary Josephson oscillations across an aperture was for a long 
time believed to be unobservable due to the strong fluctuations of the order parameter 
in the volume with dimensions of the order of coherence length. In the experiment, how- 
ever, all the difficulties have been recently overcome and a clear, unsmeared signature 
of Josephson oscillations has been found [24] . 
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Chapter 6 



Outlook 




Figure 6.1: Two examples of mesoscopic structures in which Josephson effect plays 
a role: hetero structure with a normal or ferromagnetic layer embedded between two 
superconductors: X=N,FM (a); or a quantum dot coupled to superconducting leads (b). 



The Josephson effect, predicted and discovered in the 1960-s in superconducting 
systems, opened a broad avenue of a new research area related to this phenomenon, 
which is still very active. It is impossible to mention in a short review all the interesting 
subfields which deal in one or another way with the Josephson phenomena. We give 
therefore just a few examples. 

In addition to the contributions due to Cooper pair tunneling discussed in Section 

one should take into account quasi-particle terms [2], which we did not consider. 
Those contributions are especially important for the nonstationary Josephson effect, 
i.e. in case of the finite voltage applied to the junction. Many other effects influence 
the behavior of the particle tunneling between two superconductors: various impurities, 
magnetic fields, inhomogeneities, different pairing symmetries (2} [8]. 

Many interesting phenomena arise due to the so-called proximity effect |25j: super- 
conductivity penetrates up to a certain length scale into the neighboring normal, or 
ferromagnetic material (X=N,FM as in Fig. 16.11 (a)). In the case of S-X-S heterostruc- 
ture, shown in Fig. 16.1( a) the layer between two superconductors does not need to be 
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very thin for a supercurrent to occur (for reviews see [251 EH E3 EE])- Due to the new 
experimental discoveries in superfluid Helium [23], [24], the question arises whether an 
equivalent of an S-N-S structure can be created also in these systems. 

Another interesting system is a quantum dot coupled to two superconducting leads 
shown in Fig. 16.1( b) (for review on the transport through quantum dots see for instance 
[29]). One can consider a similar arrangement in a bosonic system [30] . 

Finally, a fundamental physical problem is a nonequilibrium Josephson effect. For 
example, in a bosonic system it is easily realized, as the barrier between the two wells 
confining condensates is ramped up in a nonadiabatic way [21]. This gives rise to quasi- 
particle excitations out of the condensate [HI [32]. One can then develop a description 
in terms of the Keldysh Green's functions [32], which we mentioned in the context of 
the superconducting Josephson junction in Section El 
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